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The shape deformation of a three-dimensional axisymmetric vesicle with encapsulated filaments 
or impurities is analyzed by integrating a dissipation dynamics. This method can incorporate 
systematically the constraint of a fixed surface area and/or a fixed volume. The filament encap- 
sulated in a vesicle is assumed to take a form of a rod or a ring so as to imitate cytoskeletons. 
In both cases, results of the shape transition of the vesicle are summarized in phase diagrams 
in the phase space of the vesicular volume and a rod length or a ring radius. 

We also study the dynamics of a vesicle with impurities coupled to the membrane curvature. 
The phase separation and the associated shape deformation in the early stage of the dynam- 
ical evolution can well be explained by the linear stability analysis. Long runs of simulation 
demonstrate the nonlinear coarsening of the wavy deformation of the vesicle in the late stage. 



I. INTRODUCTION 

Amphiphilic molecules in the aqueous solvent assem- 
ble into a bilayer membrane, and form a closed vesicle. A 
typical example is an artificial lipid-bilayer, which serves 
as a model for biomembranes. Vesicles of a lipid-bilayer 
membrane have been studied as a model system of living 
cells over a couple of decades. In particular, after Helfrich 
roposed an elastic continuum model of the membrane 
many works have been made about the morphology 
of the vesicles. Even in the case of a single component, 
vesicles take various shapes, and their shapes are classi- 
fied as spherical, prolate, oblate, pear, stomatcyte and so 
on, according to Seifert et al.[^ 

Actually, biomembranes might be more complex by 
containing many components such as several different 
species of lipids or membrane proteins. Impurity lipid 
molecules and proteins can move freely within the two- 
dimensional liquid state of the membrane and perform 
diffusive motion. Via the coupling to the membrane cur- 
vature, impurities may affect the shape of the vesicle. 
Also the vesicle can enclose cytoskeletons in its interior 
space. Cytoskeletons consist of microtubles and/or in- 
termediate filaments, and form a three-dimensional net- 
work inside a cell to support the cell shape. Also for the 
cell locomotion, the elongation and the shortening of fila- 
mental proteins in the vesicle are required.]^ Thus in the 
study of living cells, the consideration of the interaction 
between lipid bilayer membranes and such encapsulated 
materials is indispensable. 

To study the role of such materials, several experiments 
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have been made. Hotani and MiyamotoQ established a 
model system where microtubules polymerize in a vesi- 
cle. A growing microtuble pushes the vesicle membrane 
at contact points and deforms it into the shape resem- 
bling a Greek letter (jj. There have been many experi- 
ments and theoretical analyses on this phenomenon. It 
was established that to realize the (/)-shape, the osmotic 
pressure inside the vesicle has to be higher than the outer 
pressure,!^, |j or the volume of the vesicle has to be fixed. 
10, H It is also found that the monopolar asymmetric (jj- 
shape has a lower energy and is stabler than the bipolar 
mirror symmetric shape. 

When the filaments in the vesicle are soft such as actin 
filaments, they bend and form a ring in the vesicle. The 
experiment of such cases was performed by Miyata and 
Hotani ||]. They observed that ring filaments in a vesicle 
push the vesicle outward and deform the spherical vesicle 
into a disk shape. 

Shape transformation caused by the phase separation 
within the membrane occurs in many important natu- 
ral phenomena such as cell locomotion, fusion, secretion, 
endocytosis, phagocytosis etc. Furthermore, various ex- 
periments have been performed on this type of shape 
changes in multi-component vesicles. One example of 
the experiments is on the transition from a biconcave 
shape of erythrocytes to a crenated one (echinocytosis) 
0, In). Another one is on a shape deformation induced 
by the phase separation of the amphiphiles composing the 
membrane p^. Also there is an experiment on a shape 
deformation caused by proteins or polymers which 
are anchored in membranes jl^. In these cases, vesicles 
transform from the spherical shape to the one with a 
budding or from a tubular shape to a pearling one. The 
effect of impurity on the vesicle morphology was first an- 
alyzed by S. LeiblerfT^. He showed that the impurities 
diffusing within the two-dimensional membrane couple 
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to the curvature and destabilize the fiat membrane. Re- 
cently, dynamical simulations of two-component vesicles 
have been performed and the coalescence process of 
the budding vesicle was obverved[^. 

In this paper, we study morphology of a three- 
dimensional axisymmetric vesicle with encapsulated fil- 
aments or embedded impurities, by solving purely dissi- 
pative dynamics. In section |ll[ the model Hamiltonian 
of a vesicle is introduced, and its dissipative dynamics 
is formulated. In section [II, the shape of the vesicle 
encapsulating a rod filament is studied. As a linear fil- 
ament elongates, vesicular shape alters from a sphere to 
a rugby-ball, and to a sphere with tubular protrusions. 
The shape transition is summarized in a phase diagram 
in the phase space of the vesicular volume and a fila- 
ment length. In section we study the case of ring 
filaments. When filaments in the vesicle are soft such as 
actin filaments, they bend and form a ring. As the fila- 
ment extends their length, the ring pushes out the vesicle 
membrane and deforms the prolate vesicle into the oblate 
one. In the case of cell division, the ring filament shrinks. 
Relations between the ring-force and the vesicle shape or 
the ring radius are studied in this case, and the morpho- 
logical transition is summarized in the phase diagram. 
In section we study the effect of impurities embed- 
ded in the membrane on the vesicle shape, and observed 
that impurities diffusing within the vesicle surface induce 
curvature instability. A linear stability analysis around 
a spherical vesicle is confirmed in our simulation of the 
vesicle under several conditions. We have further simu- 
lated the nonlinear coalescing process in the late stage. 



Whole analysis is summarized in section VI 



II. 



DISSIPATION DYNAMICS OF 
MEMBRANES 



The shape of a membrane vesicle is characterized by 
two parameters Si, {i — 1,2) such that the position vec- 
tor f of the membrane surface in three dimensional space 
is represented as r(si,S2). With this parametrization, 
the metric tensor is given by gij — diX ■ djf, and its de- 
terminant is denoted as g — det{gij}. Here di denotes the 
partial differentiation by Si as di — d/dsi. The surface 
area element is given by dA = ^d^s and the volume el- 
ement by dV = dA(f-n)/Z where n is the surface normal 
vector. 

The energy of the system consists of the bending elastic 
energy, and the surface and the volume contributions as 



E='^ j{H ^ CofdA + (tA + PV. 



(1) 



Here k is the bending rigidity, H the mean curvature 
defined as the sum of two principal curvatures Ci and 
C2 as 7J = Ci -|- C2, Co the spontaneous curvature, a 
the surface tension, and P the pressure. In a dissipative 
dynamics, the system evolves to reach minimum energy 



configuration dictated by the energy (|^). To achieve this 
purpose, we assume a purely dissipative dynamics with 
the Rayleigh dissipation function of the form 



I / \v?dA. 



(2) 



with a viscosity coefficient r?. The equation for time evo 



lution is then written as llq, 19, EQ 
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(3) 



This gives, in general, the shape evolution equation as 
1 5E 



dtf 



VnU^^Vc-i] ^aH-j] ^P)n, (4) 



where the normal component of the velocity due to the 
curvature elasticity is given by 

= - + 2AH + 4C1C2C0}. (5) 

Here A ~ g^^/^di{g^^^d^) is the Beltrami-Laplace oper- 
ator. This evolution equation (4) is explicitly derived by 
Marsili et al. jl8| in section III. A. 1-3 in their paper. Shape 
evolution generally induces variations of the surface area 
or the volume enclosed by the vesicle as[|l8[ |9[ ^ 

dA 

m 

^ = /(i^e - rr Vff - r7-ip)dA (7) 



/ H{vc - r^-^aH - r]-'^P)dA, 



(6) 



If the area of the vesicle is to be fixed, the tension a 
should be regarded as a Lagrange multiplier, and has to 
be adjusted so as to satisfy dA/dt = 0. If the volume as 
well as the surface area are to be fixed, both the tension a 
and the pressure P have to be adjusted to yield: dA/dt = 
dV/dt = 0. 

For an axi-symmetric vesicle, the position vector is con- 
veniently written in the cylindrical coordinate as (p, z), 
where the z axis is chosen to be the axis of rotational 
symmetry. As for the two parameters, we choose an ar- 
clength si = s = J ^ {dzY -\- {dpY and the angle S2 = ^■ 
The arclength is measured from the bottom of the vesicle 
on the symmetry z axis. The rotational symmetry tells 
us that the shape is independent of the second parameter 
S2 = V- Then the metric is ^Jg = p, the angle-integrated 
area element is dA = 2-Kpds and the volume element is 
dV = Tvp'^Zsds. Since the shape is independent of the ro- 
tation angle ip, we consider hereafter a section at (/? = 0, 
and the vesicle contour is given in (p, 2:)-space. The nor- 
mal and tangential vectors in this space are given as 

n^{zs,-Ps), t={ps,Zs), (8) 
and two principal curvatures are calculated to be 

Ci = Zs/p, C2 = PsZss - ZsPss- (9) 
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Here the subscript s denotes differentiation by s, for in- 
stance, Zg — dgZ. The Beltrami-Laplace operator is rep- 
resented as 

A - -d, {pd,) . (10) 

p 

The variation of the elastic energy produces only the nor- 
mal component of the velocity, because only this com- 
ponent has the physical significance. Nevertheless, we 
introduce a tangential component of the velocity, which 
corresponds to the rearrangement of the parametrization 
s, or a gauge, on a vesicle surface. In the practice of shape 
simulation, we deal with only a discrete set of points 
(called grid points, hereafter) on the membrane. Even if 
these grid points are initially prepared in equidistance, 
they will dilatate or contract locally during the shape 
evolution. In order to redistribute grid points in equal 
separation, we must add an appropriate tangential ve- 
locity Vt{s)t at a grid point of an arclength s. Since the 
total arclength L varies in time, the condition of equal 
separation of grid points are expressed by imposing the 
invariance of a relative arclength, namely d{s/L)/dt — 0. 
By this condition, we get the tangential velocity|l^, ^ 

Vt{s) = - I VnC2ds+Y / VnC2ds. (11) 

Jo ^ Jo 

The evolution of a vesicle is then obtained by numerically 
integrating the equation 

dr — 

= Vnn + Vtt. (12) 



III. ENCAPSULATION OF ROD FILAMENT 



Here we use the approximation such that around the pole 
at s = 0, the profile is almost horizontal and ds = dp. 
The integration of the force balance equation then yields 
the curvature [|| as 

H{p) ^ Ho log p/po, (15) 

where Hq and po are appropriate constants. This expres- 
sion implies that the curvature diverges at the end s = 
(p — 0) and the same type of divergence appears at the 
other end {s = L). Even though this diverging curvature, 
gives rise to only finite contribution to the bending en- 
ergy, the vesicle shape at the poles cannot be derived by 
the dynamics governed by the curvature alone. To over- 
come this difficulty, both poles are pushed outwardly un- 
til the pole-separation £ reaches the prescribed length of 
the rod filament. When the pole-separation exceeds the 
prescribed distance the poles are set free from exter- 
nal force and the distance decreases due to the curvature 
force by the nearby surface. The force applied on the 
poles are estimated from Eq. ( p^ ) by the curvature differ- 
ences between pi and p2 close to the poles as 

F = 2n.^^^M^P^. (16) 

Simulation results are normalized in the dimensionless 
form. The length is measured in unit of Rq, which is the 
radius of the sphere with a surface area of the vesicle, 
and the energy by Sttk. Then the volume is normalized 
by 47ri?Q/3, the pressure P as p = PRq/k and the force 
F as f — FRq/Attk, respectively. In the simultaion the 
spontaneous curvature Co is set to be 0. 



We consider here vesicle deformation induced by the 
elongation of an encapsulated filaments such as micro- 
tubles. The filaments are rigid and extend straight in 
the vesicle. They push and deform the incarcerating 
vesicle. There have been several studies on this prob- 
lem so far. 1^, 1^ We first reproduce known results by 
the present method. Since the encapsulated filament is 
assumed to remain straight, it exerts pressing force at 
both end points or poles at s = and s = L as FSo{s) 
and F6o{s — L), respectively. Here the invariant form of 
the delta function in the parameter space is given |l8) as 

S,(g-g)^^Al_ll, (13) 

In the equilibrium case, the forces balance equation 
around the pole is given by i]Vn{s) + FSo{s) = 0. By 
integrating over a very small area around the pole s = 
with a radius p, the second term gives simply F, whereas 
in the first term only the Beltrami-Laplace term gives 
rise to nonvanishing contribution as 

K / p'^d,{pd,H)2TTpds « 2'KnpdpH. (14) 




FIG. 1: Separation distance between two poles £ versus 
the strength of the force / for a vesicle under a pressure 
p — —80. Data obtained in the present analysis are de- 
picted by the symbol •. A solid curve represents data 
obtained by the variational method by Umeda et al. 
Results obtained by the Monte Carlo simulation ^ are 
represented by the symbol o. 

We first summerize the results of the simulations with 
a fixed osmotic pressure p. By starting from the mirror- 
symmetric shape of the vesicle with symmetric forces at 
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(c) 

FIG. 2: (a) Energy difference of the assymmetric and the 
symmetric vesicle as a function of the rod length I. The 
volume of the vesicle is fixed at u = 0.7. (b) Morphology 
phase diagram of the vesicles with various volume v and 
the rod length ^. A dashed curve indicates the cylinders 
with hemishperical caps at both ends, (c) The strength 
of the force / versus the rod length £ at a volume u = 0.7. 
The critical length for protrusion is Ic — 4.6. 

both poles, we get the vesicle with a mirror symmetry 
about the equator. The shapes as well as the length- 
force relation are shown in Figjl^, in good agreement with 
those obtained previously. ||| With a negative osmotic 
pressure, the elongation of the filament length ^ inside a 
spherical vesicle leads first to the increase of the strength 
of the force / applied to both poles, but eventually / 
reaches the maximum when two poles start to protrude. 
The maximum value of / turns out to be a little larger 
than the value obtained by Umeda et at, probably due to 
the relatively small number of grid points {N — 120) in 
our simulation. By applying asymmetric forces at both 
poles in the initial stage of the simulation, we can get an 
asymmetric shape. This is due to the fact that the asym- 
metric shape has a lower energy and is stabler than the 
mirror symmetric one, as is already discussed previously. 



One of the advantage of the present method is that 
the shape of a vesicle with a fixed volume can be stud- 
ied easily and quickly. ||] When the length of the rod is 
short, the vesicle takes a mirror symmetric shape inde- 
pendent of the initial condition. On the other hand when 
the rod is long and two protrusions appear, two distinct 
shapes, that is symmetric and asymmetric shapes, be- 
come possible. Both shapes have similar energies, but the 
energy difference between these two shapes at a volume 
V = 0.7 shown in FigJ|(a) indicates that the asymmet- 
ric shape is energetically favorable than the symmetric 
one. The morphological phase diagram is summarized in 
the phase space of the volume v and the rod length i in 
Fig.^(b). Some typical shapes are shown therein. The 
phase boundary marked by / = represents the v — £ 
relation for the free prolate shape without any force at 
poles. In order to decrease the volume under the condi- 
tion of a fixed surface area, the vesicle changes its shape 
from sphere to prolate. In fact, the v — £ relation can 
be well explained by assuming that the vesicle takes the 
shape of a cylinder with hemispherical caps at both ends. 
The radius of the sphere r and the length of the cylin- 
der £ ~ 2r should be determined from the conditions of 
constant surface area Anr^ + 2-Kr{t — 2r) = Att and the 
constant volume 47rr^/3-t-7rr'^(£ — 2r) = Airv/i. The v — £ 
relation so determined is denoted by a dashed curve in 
Fig.^(b). For a cylinder part, one of the curvature van- 
ishes, C2 = 0. In the actual situation the vesicle lowers 
its elastic energy by selecting an appropriate value of C2 , 
and thus the phase boundary f — deviates a little from 
the curve expected for the simplified cylindrical shape. 

When the force is applied on two poles of the vesicle 
with a fixed volume, the rod length £ increases. For in- 
stance at V = 0.7, the strength of the force / increases 
as £ as shown in Fig|2|(c). Initially, as the filament elon- 
gates, / increases rapidly until the poles begins to push 
out. Once one or both poles start to protrude, the incre- 
ment of / gets mild. As the rod gets long, the tubular 
protrusion becomes thin with a small diameter, and the 
central bulb gets more rounded. 

The shape of the vesicle under a strong force might 
be approximated by the form of a cylindrical tube at- 
tached to a spherical bulb. In order to fulfil the volume 
conservation, the sphere covers the volume 47rv/3 such 
that its radius is just v-^^^. Since the sphere gives only 
the area Attv"^^^, the remaining area is covered by the 
fine cylindrical tube with an extension £' and the radius 
r « 2(1— For £' — > 00, the volume of the cylinder 
vanishes irr'^f ^ 1/f ^ 0. Therefore, the rhs region of 
the phase diagram in Fig. 2(b) extends to infinity except 
at t; = 1. 



IV. ENCAPSULATION OF RING FILAMENT 

When the filaments encapsulated are soft like actin 
filaments, they tend to bend and to form a ring in a 
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vesicle. As their degree of polymerization or their length 
increases, the ring radius expands. On the contrary, in 
the case of cell division, an actin filament ring shrinks. 

The force exerted by the ring at z = can be estimated 
by assuming the force balance: riVn{z = 0)+F/2TTrS(z) — 
0. Here F is the strength of the total force exerted by 
the ring of radius r at z = 0. Similar to Eq.(14), an area 
integration around z = can be performed with p fa r 
and ds — dz. Then it yields the curvature as 

H{z) = Ho- -^z. (17) 

where Hq is the curvature at z = 0. In this case, Hq 
remains finite. The total force applied by the ring at 
z = is estimated from Eq.(p^ by taking the curvature 
difference as 



F ^ 2'KrK 



Ho - H{zi) 



z\ 



(18) 



We now simulate the shape variation of the vesicle in- 
duced by the force exerted by the ring filaments. Every 
physical quantity is scaled in the dimensionless form as 
before. For a given volume v = 0.7, the free vesicle has a 
prolate shape with a radius about r « 0.42 at the equator 
z = 0, and the other curvature C2 is negative. As the ring 
radius increases, the bending energy increases, as shown 
in FigJ^(a). Correspondingly, the force calculated by the 
formula Eq.(18) also increases, as shown in Fig.||(b). The 
dotted curve shows the result obtained by numerically 
differentiating the energy shown in Fig.^a) by the ra- 
dius r. This indirect calculation of the force agrees well 
with the force obtained directly from the simulation until 
the force passes the maximum for the prolate shape and 
the minimum for the oblate shape. 

With a small increment of the radius, the curvature C2 
at the equator is still negative. On increasing the ring 
radius, the central zone expands to have two positive 
curvatures, still keeping the two poles far apart as in 
the prolate shape. As the central radius increases, the 
force to expand the radius starts to weaken. At a large 
enough radius as r > 1.14, the prolate becomes unstable 
and the vesicle jumps to an oblate shape. The oblate 
at r = 1.23 is, in fact, a shape with a local minimum 
energy without a force. It has a bending energy being 
a little more than that of the prolate shape, and thus is 
a metastable shape. On the further increase of the ring 
radius, the shape changes to that of a flying saucer with 
a large energy cost. On decreasing the ring radius from 
that of the free oblate shape, the dips at poles deepen. At 
a small ring radius, another dips develope at the equator, 
too. Actually, the bending energies of the prolate and the 
oblate shapes cross at a critical radius Vc = 0.96 for the 
fixed volume v = 0.7, and the first-order shape transition 
should take place there. The morphology phase diagram 
with prolate and oblate shapes is shown in the phase 
space of the volume v and the ring radius r as shown in 
Fig.|(c). 
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FIG. 3: (a) Bending energy E versus ring radius r of 
a vesicle with a volume v = 0.7 for the prolate (filled 
circles) and the oblate (filled triangles) shapes. Lines 
are guides for the eyes, (b) The strength of the radial 
force / versus the radius r at the equator. A dashed 
lines indicate the strength of the force induced by the 
derivative of the bending energy, (c) Morphological phase 
diagram in v and r space. A dashed curve indicateds the 
upper and lower limiting shapes (see text). The critical 
radius is Tc = 0.96. 



When the ring shrinks in the middle of the prolate 
shape from the natural radius r = 0.42, as in the case of 
cell division, the two parts, the northern and the southern 
area, get round, as shown in Fig. 3(a) and (b). Since the 
surface area and the volume of the vesicle are fixed, the 
central radius can vanish only for a vesicle with a normal- 
ized volume less than l/^/2 m 0.71. When the volume is 
larger than this value, the shrinkage stops intermediately 
with a shape of two connected hemispheres. 

When the volume is smaller than this value, the shape 
of the vesicle at r = is equivalent to the two connected 
equilibrium shapes with half of the volume. As shown in 
Fig.^(a), when the volume of the vesicle is 0.7, the vesi- 
cle approaches to the shape with two connected almost 
spherical shapes with v = 0.99, where the bending energy 
of the vesicle is about 2. 

When the ring shrinks or expands in a vesicle with a 
fixed volume and a fixed surface area, there is a certain 
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limit for the ring radius r determined solely by the geo- 
metric reason. For a ring with a minimum radius r, the 
vesicle takes the symmetrical semi-spherical shape; the 
radii of upper and lower spheres are the same as R and 
the separation between two centers is 2d. For a ring with 
a maximum radius r, the vesicle again takes the symmet- 
ric semi-spherical shape with the separation between two 
centers equal 2c? being negative. The area conservation 
gives the condition AttR{R + d) = An, and the volume 
conservation 2tt{R + df{2R - rf)/3 = A-nv/i. The ra- 
dius at the equator is given by r = V2 — i?^^. When 
d — R — l/\/2, the vesicle consists of two complete 
spheres with r — and the volume is equal v = 1/^/2. 
When d = —R = —oo, the vesicle becomes flat with a 
vanishing volume w = 0, but the radius at the equator 
remains finite r = V2. The upper and lower limits of 
r for arbitrary volume v can be easily obtained numeri- 
cally, and are shown by a dashed curve in Fig.^(c). At 
V = 0.7, the upper limit of the radius r thus obtained 
agrees well with the result of dynamical simulation. 



V. IMPURITY DIFFUSION 

There are usually various ions or proteins dispersed 
in a membrane, and they affect the shape of vesicles. 
When the impurities are homogeneously distributed, the 
vesicle takes a shape determined by the curvature elastic- 
ity with some modification in the spontaneous curvature. 
But when the impurities are coupled to the local curva- 
ture of the membrane, the homogeneous distribution may 
become unstable. The phase separation and the shape in- 
stability are induced simultaneously. We here study this 
phenomenon analytically and numerically. 



A Dynamics driven by impurities 

Since the homogeneously distributed impurities only 
modify the spontaneous curvature and the surface ten- 
sion cr, we set the average concentration to be zero. The 
local concentration fluctuation of intercalated molecules 
from the average is denoted by ^{f,t). The interac- 
tion of intercalated molecules with the phospholipidic 
constituents of the membrane might be summarized in 
the form, 



where the homogeneous part is expanded up to the fourth 
order as 



(21) 



Since the average fluctuation of impurity concentration is 
set to zero, there is no linear term of in the free energy 
/(<i>). The evolution dynamics of the membrane, Eq.(|^) 
now contains terms due to the impurity as 

vn ^vc + 77-^A{2CiC2$ - A$} - ,y-i{|(a,$)2(Ci - C2) 
+ i^'i>^ + ^^'^)H}-Tj-^aH -7]-^?. (22) 

The conservation of the total number of impurity 
molecules on the membrane leads to the restriction 
Jd^Sy/g^ = 0. Therefore, the time-evolution of the 
impurity concentration fluctuation $ should satisfy 



«„). (23) 



This is satisfied by describing the evolution in terms of 
the diffusion equation as 



(24) 



In a simulation where a membrane is discretized into grid 
points, the supporting grid points have a tangential ve- 
locity in addition to the physical normal velocity, so as 
to keep the grid separation equidistant. Accordingly, the 
concentration fluctuation field $ should include this con- 
tribution given as 

^ d^m) ^^^dr ^ ^^^^ 

dt dt ^ ' 

Therefore, the flnal evolution for <i> is described as 

dt^ + <S>Hv„-vtds^ = DA{-BA^ + A2^ + Ai'A>^-AH). 

(26) 

We introduce the following dimensionless parameters 
hereafter. 



0.2 



A2 

nRl' 



a4 



A 

kRc 



, d — tiDRq, Co = CqRo- 



(27) 



-A 



d^s 



(19) 



with a coupling constant A. The free energy of the impu- 
rity molecules is assumed to be Landau-Ginzburg form 



F 



d^s^[-B{V^)' 



fm 



(20) 



B Linear stability analysis 

In this section we study the stability of the spherical 
vesicle with homogeneously distributed impurities under 
a small shape fluctuation and the impurity inhomogene- 
ity. For simplicity, we consider the case with a fixed 
pressure P, and thus the volume of the vesicle can be 
altered. In the unperturbed situation, the vesicle has a 
radius Rq. Assuming the small deformation of the vesi- 
cle and no overhang configuration, the vesicular position 
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vector r and the concentration of impurities in dimen- 
sionless form are single-valued functions of the polar 
and the azimuthal angle (p in the spherical coordinate. 
These small deviations are expanded in terms of spheri- 
cal harmonic function y^, and we note that each mode 
is independent in the linear stability analysis. Therefore, 
one considers only the Z-pole mode, and the membrane 
position vector r and the impurity concentration <j) in 
dimensionless form are written as 



(29) 



where i5 and C, describe the time-dependent coefficients 
of fluctuations and is the unit radial vector. Since 
the average concentration is set to zero, there is no un- 
perturbed term in the expansion of the concentration (f>. 
From Eqs.(p2|) and (^), one obtains the linear dynamics 
of the perturbation 6^, (1^ as 



r{e, v,t) = {l + Sl{t)Yi{0, p)}er{0, (28) 



dt 
dt 



(/ + 2)il - 1) |-?(/ + 1) + + coj Si + {2A + Mil + 1)} CL 
d{-l{l + 1) [Mil + 1) + 02] CL + + + - l)Sl} . 



(30) 
(31) 



The conservation of the surface area determines the surface tension a. With a small shape deformation 6, the variation 
of a is easily shown to be quadratic in 5, and it does not affect the linear stability analysis. Eqs. ( ^ ) and ( ^l| ) can 
be reduced to matrix from as follows. 



-{1 + 2){l - 1) {1(1 + 1) - f - Co} {2 + l{l + 1)}A 

dl{l + + 2){l - 1)A -dl{l + l){bl(l + 1) + 02} 



Sm 



(32) 



Note that the time-development matrix M is independent of m. The eigenvalues /3+,/9_ and eigenvectors {S±,(±) of 
the matrix M define the time evolution of the deviations Cm as follows. 

( ) = exp(/3+i) ^+ ^ + A_ cxp(/3_t) (33) 

/3± = ^{Mrr + A/00 ± ^(Af„ - M00)2 + 4Mr4>M4,r}, 



c± 



In Fig.Q we show the relationship between the growth 
rate f3+,P- and mode number I for both cases without 
A = (solid curves) and with the coupling A ^ (dashed 
curves). 

With a negative 02 , as in Fig.^(b) impurities can phase 
separate spontaneously, i.e. even if there is no coupling to 
the curvature of the membrane, M^^ > for Z* > 1 when 
2b < —02. Then the most unstable mode Za=o inte- 
ger close to l/2(-l-t-^l - 202/6). When2& < -02 < 6b, 
the most unstable mode is /* = 1. The impurities tend 
to phase separate into two domains. With a small cou- 
pling to the membrane curvature (A 7^ 0), this mode 
induces the translation of the vesicle without deforma- 
tion. When 6b < —02, some modes with I* > 2 becomes 
unstable. With a small coupling A, the vesicle deforms 
accordingly. 

To confirm the linear stability analysis and further to 
elucidate the nonlinear effect, we simulate the dynamical 
evolution of the vesicular shape and the impurity diffu- 
sion. We assume axisymmetry in the simulation, and 
thus only m = modes are relevant here. The parame- 
ters chosen are b = 0.1, 02 = —1, 04 = 1, A = 3.0, d = 



1, Co = 0. The pressure p is fixed to be p = 0. The most 
unstable mode is expected to be Z* = 6. 

The initial shape is spherical, and the impurities are 
distributed randomly whose Fourier amplitudes up to 
32nd mode Co for ^ < 32 have uniformly distributed 
between —0.001 and 0.001. In the course of simulation 
the vesicle shape developes the wavy pattern with 6 do- 
mains of impurities, as shown in Fig.^(a). Impurities 
accumulate to the convex portion of the vesicle due to 
the positive coupling A > 0. This number corresponds 
to the most unstable mode with I* — 6. As the time 
evolves, these domains coalesce each other. In last stage 
of Fig. 5(c) the number of the domains becomes 4. 

With a positive 02, accumulation of impurities costs 
energy and they try to diffuse away to be homogeneous 
if there is no coupling to the curvature of the membrane. 
However, due to the coupling between the concentration 
variation and the membrane curvature, the homogeneous 
distribution can be unstable and the phase separation 
is possible for sufhciently strong coupling as shown in 
Fig.^(a) . The linear stability analysis shows that a spher- 
ical vesicle becomes unstable if the coupling is stronger 
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(a) 



(b) 




(a) 



(b) 



(c) 



FIG. 4: Growth rate /3 as a function of mode number 
I in the hnear stability analysis for P=0, Co = 0, (a) 
for 02 > 0, (b) for a2 < 0. Solid curves represent the 
growth rates for A — 0, (3+ (filled circles) and /3_ (filled 
triangles). Dashed curves represent the growth rates for 
A 7^ 0, /3+ (filled squares) and f3- (filled diamonds). The 
parameters are b — 0.1, 02 = 1 or —2, 04 = 1, A = 
2.0 d= 1, Co = 0,p = 0. 




■3.0 3.0 



FIG. 6: Time evolution of the vesicle with impurities, 
the parameters are 5 = 0.1, 02 = 1, = I, X — 3.0 d — 
1, Cq — 0,p = 0. The most unstable mode is Z* = 6 




2.2 2.2 



(a) 



(b) 



(c) 



FIG. 7: Time evolution of the vesicle with impurities, 
the parameters are b = 0.1, 02 = 1, 04 = 1, A = 2.5 d = 
1, Co = 0,_p = 0. The most unstable mode is P = 5 



(a) 



(b) 



(c) 



FIG. 5: Time evolution of the vesicle with impurities, 
the parameters are 6 = 0.1, 02 = — 1, 04 = 1, A = 
3.0 d = 1, cq = 0,p = 0. The most unstable mode is 
I* = Q The time is measured in unit of tiRq/k. 



than the critical value |A*| = \/18b + 3a2/2. For the pa- 
rameter choice b — 0.1, a2 — 1, a4 = 1, d = 1, co = 0, 
the critical coupling takes the value A* ~ 1.09. When the 
coupling is stronger than A* such as A = 3.0, the most 
unstable mode has /* = 6, in good agreement to the early 
stage of the simulation result as shown in Fig||(a) . 

By decreasing the coupling such as A = 2.5, the number 
of the most unstable mode decreases to P = 5, in agree- 
ment to the simulation result shown in Fig.|^(a). When 
the correlation length b gets smaller such as 6 = 0.05, the 
number of the most unstable mode increases as I* ~ 9. 
In the simulation the vesicle becomes wavy with 8 do- 
mains in the early stage as shown in Fig.||(a). The result 
has the same tendency with the linear stability analysis. 
In all the cases studied, the coarsening takes place as the 
time elapses, and the number of domains decreases. 

We can study the effects of the pressure p or the spon- 
taneous curvature co by linear analysis. According to 
Eg. (p2[) , the positive p or co increase the number of 
the most unstable mode /*. On the contrary the neg- 
ative p or Co decrease I*. For example, we simulated 




(a) 



r ' 




I- 


Q.S^ 1 

.14 

(b) 




-2.6 2.6 



(c) 



FIG. 8: Time evolution of the vesicle with impurities, 
the parameters are b — 0.05, 02 = 1, 04 = 1, A = 3.0 d = 
1; Co 0,p = 0. The most unstable mode is I* —9 








t=0 
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28 



■2.5 2.5 



(a) 



(b) 



(c) 



FIG. 9: Time evolution of the vesicle with impurities, 
the parameters are b = 0.1, 02 = 1, 04 = 1, A = 3.0 d = 
1, Co = 0,p = —40. The most unstable mode is Z* = 5 
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the shape of the vesicle with the parameters such as 
p = -40, b = 0.1, 02 = 1, a4 = 1, A = 3.0 d=l, cq = 
as shown in Fig.^. Then the linear analysis gives the 
most unstable mode I* — 5, in fair agreement with the 
simulation. From Fig. 9(b) and (c), the coarsening does 
not seem to take place under the pressure. This might 
be related to the resistence to the volume change at a 
negative p. 

In almost all cases the results of the simulation in the 
early stage agree with the linear stability analysis. More- 
over, the phase separation domains coalesce each other, 
as the time evolves. This nonlinear coalescing process is 
also observed in the previous study p^. 



VI. CONCLUSION 

We performed simulations of the axisymmetric vesi- 
cle shape change caused by the length change of an en- 
capsulated ring or rod filaments. The constant volume 
simulation of the vesicle with a rod encapsulation shows 
that the (j) shape consisting of a spherical and a tubu- 
lar portions is attributed to the restriction of the volume 
and the surface area conservations. The same restriction 
limits the extention of the ring radius in the case of the 
ring-filament encapsulation. The shape transitions of the 
vesicle are summarized in phase diagrams. Fig. 2(b) and 
3(c), in the phase space of the vesicular volume and a 
rod-filament length or a ring- filament radius. ^, ^ |j. 

We also performed simulations on the dynamics of 
the vesicle shape caused by the phase separation of the 
impurity concentration. Impurities which show sponta- 
neous phase separation induces periodic deformation in 
the vesicle shape. Even those impurities which do not 
give rise to spontaneous phase separation are shown to 
induce periodic deformations through the strong coupling 
of the impurity concentration to the membrane curva- 
ture. The early stage is well explained by the linear 
stability analysis, such as the dominant periodicity of 
the shape deformation. In late stage where the nonlin- 
ear effect becomes dominant, the coarsening of the shape 
periodicity takes place. Similar coasening process was 
observed in the previous study[0. The study on the 
strongly deformed vesicle due to the impurity phase sep- 



aration is now under way. 
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